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Abstract 



We report a measurement of the W boson mass based on an integrated 
luminosity of 82 pb~^ from pp collisions at ^ = 1.8 TeV recorded in 1994- 
1995 by the D0 detector at the Fermilab Tevatron. We identify W bosons by 
their decays to eu, where the electron is detected in the forward calorimeters. 
We extract the mass by fitting the transverse mass and the electron and 
neutrino transverse momentum spectra of 11,089 W boson candidates. We 
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measure M]^ = 80.691 it 0.227 GeV. By combining this measurement with our 
previously pubhshed central calorimeter results from data taken in 1992-1993 
and 1994-1995, we obtain Mw = 80.482 ± 0.091 GeV. 



Typeset using REVT^]X 



5 

In the standard model of the electroweak interactions (SM), the mass of the W boson is 
predicted to be 

^ V V2Gf ) sin^vi/Vl- Ar 

In the "on-shell" scheme [|I| cosOw = Mw/Mz, where Mz is the Z boson mass. A mea- 
surement of Mw, together with Mz 0, the Fermi constant (Gp), and the electromagnetic 
coupling constant (a), experimentally determines the weak radiative corrections Ar. Com- 
pared to the formulation in [jl| where a was defined at = 0, we have absorbed purely 
electromagnetic corrections into the value of a by evaluating it at = M|. The dominant 
SM contributions to Ar arise from loop diagrams involving the top quark and the Higgs 
boson. If additional particles coupling to the W boson exist, they also give contributions to 
Ar. Therefore, a measurement of M^/ is a stringent experimental test of SM predictions. 
Within the SM, measurements of Mw and the mass of the top quark constrain the mass of 
the Higgs boson. 

We report a new measurement of the W boson mass using electrons detected at forward 
angles. We use 82 pb~^ of data recorded with the D0 detector during the 1994-1995 run 
of the Fermilab Tevatron pp collider. This forward electron measurement, in addition to 
increasing the statistical precision, complements our previous measurement with central 
electrons because the more complete combined rapidity coverage gives useful constraints 
on model parameters that reduce the systematic error. A more complete account of this 
measurement can be found in Ref. 

At the Tevatron, W bosons are produced through qq' annihilation. W ^ eu decays are 
characterized by an electron with large transverse energy {Et) and significant transverse 
momentum imbalance {j^t) due to the undetected neutrino. The particles recoiling against 
the W boson are referred to collectively as the "underlying event." 

The D0 detector |Q consists of three major subsystems: a tracking detector, a calorime- 
ter, and a muon spectrometer. The tracking detector consists of a vertex drift chamber, a 
central drift chamber (CDC), and two forward drift chambers (FDC). The CDC covers the 
pseudorapidity [rj = — In (^tan |^ where 6 is the polar angle) region |?7| < 1.0. The FDC 
extend the coverage to \r]\ < 3.0. The central calorimeter (CC) and two end calorimeters 
(EC) provide almost uniform coverage for particles with \r]\ < 4. 

At the trigger level, we require > 15 GeV and an energy cluster in the electromagnetic 
(EM) calorimeter with > 20 GeV. The cluster must be isolated and have a shape 
consistent with that of an electron shower. 

During event reconstruction, electrons are identified as energy clusters in the EM 
calorimeter, which satisfy isolation and shower shape cuts and have a drift chamber track 
pointing to the cluster centroid. We determine forward electron energies by adding the 
energy depositions in the calorimeter within a cone of radius 20 cm, centered on the clus- 
ter centroid. The electron momentum (p(e)) is determined by combining its energy with 
the direction obtained from the shower centroid position and the drift chamber track. The 
trajectory of the electron defines the position of the event vertex along the beamline. 

We measure the sum of the transverse momenta of all particles recoiling against the W 
boson, ut = J2i Ei sin Oiitf, where Ei is the energy deposition in calorimeter cell i, uf is the 
unit transverse vector pointing from the beamline to the cell center, and 6i is the polar angle 
defined by the cell center and the event vertex. The ut calculation excludes cells occupied 
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by the electron. The transverse momenta of the neutrino, Pt{i^) = ~Pt{g) ~ ut, and the W 
boson, Pt{W) = —Uti are inferred from momentum conservation. 

We select a W boson sample of 11,089 events by requiring pt{v) > 30 GeV, ut < 15 GeV, 
and an electron candidate with 1.5 < |?7| < 2.5 and Pt{g) > 30 GeV. 

We extract the W boson mass from the spectra of the electron Pt{g), neutrino Pt{i^), 



and the transverse mass, = y 2pT(e)pT(i^)(l — cos A0), where A(f) is the azimuthal 
separation between the two leptons. For each spectrum we perform a maximum likelihood 
fit to the data using probability density functions from a Monte Carlo program. We model 
the production dynamics of W bosons and the detector response to predict the spectra. The 
rriT, Pt{€.), and Pt{^) spectra have quite different sensitivities to the W boson production 
dynamics and the recoil momentum measurement. By performing the measurement using 
all three spectra we provide a powerful cross-check with complementary systematics. 

Z bosons decaying to electrons provide an important control sample. We use them 
to calibrate the detector response to the underlying event and electrons, and to constrain 
the model for vector boson production used in the Monte Carlo simulations. We trigger 
on Z ^ ee events having at least two EM clusters with Et > 20 GeV. We accept Z — > ee 
decays with at least one forward electron with 1.5 < |?7| < 2.5 and another forward or central 
(|?7| < 1.0) electron. A central electron is required to have pr > 25 GeV but is allowed not 
to have a matching drift chamber track. The forward electron candidate is required to have 
Pt > 30 GeV and a matching drift chamber track. This selection accepts 1,687 Z boson 
events. 

We use a fast Monte Carlo program developed for the central electron analyses 0,0, 
with modifications in the simulation of forward electron events. The program generates 
W and Z bosons with the i] and pr spectra given by a calculation ^ which used soft 
gluon resummation and the MRST [|| parton distribution functions. We use the relativistic 
Breit-Wigner line shape with mass-dependent width, skewed by the mass dependence of the 
parton luminosity. The measured W and Z boson intrinsic widths are used. The angular 
distribution of the decay electrons includes a pT(W^)-dependent C(as^) correction |0. The 



program also generates W —>■ ei"y [^, Z — > ee7 |Tl]], and W ^ tv ^ evVv decays. 

The program smears the generated p[e) and ut vectors using a parameterized detec- 
tor response model and applies inefficiencies introduced by the trigger and offline selection 
requirements. Backgrounds are added to the Monte Carlo samples. The parameters are 
adjusted to match the data. 

The electron energy resolution {5E/ E) is parameterized by calorimeter sampling, noise, 
and constant terms. In the Monte Carlo simulation of forward electrons we use a sam- 



pling term of 15.7%/v/£'/GeV, derived from beam tests [|12|. The noise term is determined 



by pedestal distributions taken from the W boson data. We constrain the constant term 
to cec = l-Oli;o% by requiring that the predicted width of the dielectron invariant mass 
spectrum be consistent with the Z boson data. 

Beam tests show that the electron energy response of the end calorimeter can be parame- 
terized by a scale factor a-^c and an offset 5ec- We determine these in situ using Z ee de- 
cays [g. For forward electrons we obtain 5ec = -O.liO.7 GeV and oec = 0.95179±0.00187 
by fitting the observed mass spectra while constraining the resonance masses to the Z bo- 
son mass. The uncertainty on uec is dominated by the finite size of the Z boson sample. 
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Figure |l| shows the observed mass spectra from the dielectron samples and the hne shapes 
predicted by the Monte Carlo simulation for the fitted values of cec, «eC) and Sec- The 
background was determined from a sample of events with two EM clusters failing electron 
quality cuts. The calibration of the electron polar angle uses muons from pp collisions 
and cosmic rays to calibrate the drift chambers, and Z ee decays to align the EC with 
the drift chambers. 




60 80 100 120 60 80 100 120 



m(ee) (GeV) m(ee) (GeV) 

FIG. 1. The dielectron invariant mass distribution of the CC/EC (left, x^/dof = 14/19) and 
EC/EC (right, x'^/dof = 12/17) Z boson data (•). The solid line shows the fitted signal plus 
background shape and the small hatched area shows the background. The fitting window is 
70 < m(ee) < 110 GeV. 



We calibrate the response of the detector to the underlying event using the Z boson 
data sample. In Z ^ ee events, momentum conservation requires priee) = —ut, where 
Priee) is the sum of the two electron px vectors 0]. We constrain the detector response i?rec 
using the mean value of the priee) + ut projection on the inner bisector of the two electron 
directions. Z boson events with two forward electrons give a recoil response measurement 
that's consistent with the measurement performed in the central dielectron analysis 

The recoil momentum resolution has two components: a stochastic term, which we model 
as Srec/ \/pT/GeV ; and the detector noise and pile-up, which we model using the scaled ^/r 
from random pp interactions . We constrain the model by comparing the observed rms of 
Priee) + ut/ Rrec with Monte Carlo predictions. The model tuned for the central electron 
analysis gives a good description of the PT(ee) + Ux/Rrec distributions for our Z boson 
event sample. Figure ^ shows the comparison between the W boson Monte Carlo and the 
data of the projection of recoil momentum on the direction of the forward electron {u\\) and 
on a direction perpendicular to the electron momentum {u±). 



>.0.2F 




U|| (GeV)' ' (GeV)^^ 

FIG. 2. Probability distributions of (left, xVdof = 25/15) and m_l (right, xV^of = 14/15) 
for the forward W boson data (•) and the Monte Carlo simulation ( — ). 



Backgrounds in the W boson sample are due to W tu evvv decays (1%, included 
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in the Monte Carlo simulation), hadrons misidentified as electrons (3.64±0.78 %, determined 
from the data) and Z — ee decays (0.26±0.02 %, determined from herwig Q and geant 



1^ simulations). Their shapes are included in the probability density functions used in the 



fits. The results of the fits to the mr, Prie), and Pt{i^) distributions are shown in Fig. y 
and Table I. 

We estimate the systematic uncertainties in (Table II) by varying the Monte Carlo 
parameters within their uncertainties. We assign an uncertainty that characterizes the range 
of variations in Mw obtained when employing several recent parton distribution functions: 
MRST, MRS (A') 0, MRSR2 |T§, CTEQ3M [|T71, CTEQ4M |TB[ and CTEQ5M [|T|. We 
have checked that the pdf's reproduce the ri{e) distribution for the W bosons well [^. We 
allow the Pt{W) spectrum to vary within constraints derived from the priee) spectrum of 
the Z boson data 0] and from Aqcd • Smaller uncertainties in Mw are due to the removal 
of the cells occupied by the electron from the computation of ut, and the modeling of trigger 
and selection biases H]. The uncertainty due to radiative decays contains an estimate of the 



effect of neglecting double photon emission in the Monte Carlo simulation ||20|| . 

The total systematic errors are shown in Table I. The good agreement of the three 
fits shows that our simulation models the W boson production dynamics and the detector 
response well. Fits to the data in bins of luminosity, 0(e), 77(e), and ut and with changes 
to the fit window show no evidence of systematic biases. 
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FIG. 3. Spectra of 
(a) rriT, (b) pri^), and 
(c) Pt{^) from the 
data (•), the fit ( — ), and 
the backgrounds (shaded). 
The arrows indicate the fit 
windows. 
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As a consistency check, we fit the transverse mass distribution of the Z — >■ ee events. 
We retain one electron in the EC and ignore the energy of the other electron (in the CC 
or EC). The fitted Z boson mass (Fig. g) is 92.004 ± 0.895 (stat) GeV for the CC/EC 
sample, and 91.074 ± 0.299 (stat) GeV for the EC/EC sample. The combined mass is 
91.167 ± 0.284 (stat) GeV. These results are consistent with the input Z boson mass we 
used to calibrate the detector response. 

We combine the mr, PT(e), and Pt{v) measurements of My/ using a full covariance 
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TABLE I. The fitted values and errors of tfie forward W boson mass measurements in GeV. 
The confidence level (C.L.) is given by the probability of the fit. 



Fit 


Mass 


Stat. 


Syst. 


Total Error 


C.L. 


rriT fit 


80.757 


0.107 


0.204 


0.230 


81% 


pxie) fit 


80.547 


0.128 


0.203 


0.240 


8% 


pxii^) fit 


80.740 


0.159 


0.310 


0.348 


33% 



matrix that takes into account correlations between all the parameters describing the W 
boson production model and detector response, as well as the statistical correlations. The 
combination of all three forward electron measurements yields a W boson mass of = 
80.691 ± 0.227 GeV. We also combine the three central electron measurements with 
the three forward W boson mass measurements to obtain the combined 1994-1995 data 
measurement of Mw = 80.498 ± 0.095 GeV. The is 5.1/5 dof, with a probabihty of 41%. 
Further combining this with the measurement from the 1992-93 |^ data gives the 1992-95 
data measurement of = 80.482 ±0.091 GeV. This measurement subsumes all previously 
published measurements of the W boson mass by D0. 




^ m^(e,e^v) (GeV) m^(e,e^v) (GeV) 

FIG. 4. Spectra of the Z boson transverse mass from the CC/EC data (left) and the EC/EC 
data (right). The superimposed curves show the maximum likelihood fits and the hatched regions 
show the estimated backgrounds. 



From Eq. 0, using a(M|) = (128.88 ±0.09)"^ we find Ar = -0.0322 ±0.0059, which 
establishes the existence of loop corrections to M]^ at the level of five standard deviations. 
Taken together with our measured top quark mass {rrit = 172.1 ± 7.1 GeV our value of 
the W boson mass is consistent with measurements by CDF p3[ and the LEP experiments 
and with the SM prediction for a low mass Higgs boson (i.e. rriH < 100 GeV), and is in 
even better agreement with predictions |2^ in the MSSM framework. 

We thank the Fermilab and collaborating institution staffs for contributions to this work, 
and acknowledge support from the Department of Energy and National Science Foundation 
(USA), Commissariat a L'Energie Atomique (France), Ministry for Science and Technol- 
ogy and Ministry for Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments 
of Atomic Energy and Science and Education (India), Colciencias (Colombia), CONACyT 
(Mexico), Ministry of Education and KOSEF (Korea), and CONICET and UBACyT (Ar- 
gentina). 
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TABLE II. Uncertainties in the combined mj-, PT(e), and pt(j^) W boson mass measurement 
in MeV, for the forward sample (first column), and the combined central and forward 1994-1995 
sample (second column). 



Source 


Forward 


Forward + Central 


W boson statistics 


108 


61 


Z boson statistics 


181 


59 


Calorimeter linearity 


52 


25 


Calorimeter uniformity -- 8 


Electron resolution 


42 


19 


Electron angle calibration 


20 


10 


Recoil response 


17 


25 


Recoil resolution 


42 


25 


Electron removal 


4 


12 


Trigger and selection bias 


5 


3 


Backgrounds 


20 


9 


Parton distribution functions 


17 


7 


Parton luminosity 


2 


4 
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